Experimental Observation of Quantum Correlations in Modular Variables 
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We experimentally detect entanglement in modular position and momentum variables of photon 
pairs which have passed through D-slit apertures. We first employ an entanglement criteria recently 
proposed in [Phys. Rev. Lett. 106, 210501 (2011)], using variances of the modular variables. We 
then propose an entanglement witness for modular variables based on the Shannon entropy, and 
test it experimentally. Finally, we derive criteria for Einstein-Podolsky-Rosen-Steering correlations 
using variances and entropy functions. In both cases, the entropic criteria are more successful at 
identifying quantum correlations in our data. 



I. INTRODUCTION 

Quantum entanglement manifests itself in both dis- 
crete and continuous variable quantum systems. Identi- 
fying and characterizing entanglement in continuous vari- 
able systems is quite challenging, due to the large Hilbert 
space structure and variety of measurements that can be 
performed. This is especially true if the quantum state 
of interest is not Gaussian [J Q. Recently, an entan- 
glement criteria was introduced for modular variables. 
These correspond to the integer and remainder compo- 
nents of a continuous variable, such as the position x or 
momentum p Q. Modular variables naturally arise in 
interference effects [|J , and provide a natural way to con- 
vert a continuous variable into a discrete component and 
a range-limited continuous component. As such, they 
should find application in continuous-variable quantum 
information processing. 

Entanglement witnesses using modular variables, such 
as that introduced in Ref. [3[, can identify entangle- 
ment in interesting quantum states. They are particu- 
larly well suited for bipartite states that describe non- 
local interference of two-particle wave packets. A large 
series of experiments studying exactly this kind of non- 
local interference using entangledphoton pairs has been 
performed over the last decade [5 8]. Using the spatial 
correlations produced by spontaneous parametric down- 
conversion (SPDC) H, it is possible to engineer a va- 
riety of two-photon entangled states by controlling the 
pump beam 0, H, [Iol - [l2j . In particular, if each photon 
passes through a D-slit aperture, the spatial correlations 
inherent in the SPDC process can be used to produce 
entangled Z?-dimensional qudits Q. We also note that 
the spatial entanglement present in two-photon systems 
is also intricately related to quantum imaging and ghost 
diffraction |l,[l|-[l6|. 

Higher-dimensional quantum systems are interesting 
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for quantum information tasks such as quantum cryptog- 
raphy [It], Ell an d bit commitment and several stud- 
ies have been performed for spatially entangled photons 
[2(il - |24| . Qudits also present novel possibilities for funda- 
mental tests of quantum theory [25rf27| , and to maximize 
the secure information capacity of entan gled photons in 
the presence of experimental limitations [28| . 

Though there have been a number of experiments in- 
vestigating aspects of entanglement in spatial qudits, use- 
ful entanglement witnesses have been lacking. Previ- 
ous experiments have used the conditional nature of the 
fringe pattern, which depends upon the spatial coordi- 
nates of both photons and only appears in the two-photon 
coincidence counts, to substantiate the non-separability 
of the system 0, [2^, [3(|. Discrete methods, based on 
the Schmidt decomposition, have been used to experi- 
mentally infer the entanglement in this system under the 
assumption that the state is pure fl0| . 

Here we produce pairs of photons in spatially entan- 
gled states of dimension D x D and use witnesses based 
on modular position and momentum variables to iden- 
tify quantum correlations. We first test an entanglement 
witness similar to that derived in Ref. 0] to experi- 
mentally identify modular entanglement in the near-field 
(x) and far-field (p) variables, and obtain legitimate vi- 
olation only for D = 2. Then, motivated by previous 
work in continuous variables [3l| , we derive an entangle- 
ment witness using the Shannon entropy function, and 
observe that it is more successful at identifying quan- 
tum correlations in our experimental results. We fur- 
ther ask if the quantum correlations are sufficient to 
violate Einstein- Podolsky- Rosen-steering (EPR-steering) 
criteria for modular variables. EPR-steering is a quan- 
tum correlation that is strictly stronger than entangle- 
ment [H, HH , and can be a resource for one-sided device 
independent quantum key distribution |34| . In section 
IVl we use the results of Q and our previous results to 
obtain EPR-steering inequalities based on variances and 
the Shannon entropy. The entropic EPR-steering wit- 
ness identifies entanglement in our data. An advantage 
of these witnesses is that they are valid for any initial 
bipartite quantum state. 
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FIG. 1: (Color online). Experimental setup, (a) State prepa- 
ration. A CW laser generates collinear photon pairs in a BiBO 
crystal. Spherical lenses SL1, SL2 and cylindrical lens CL are 
used to increase the efficiency of SPDC and image the spatial 
correlations of the photons at the D-slit (D = 2, 3, 4) plane, 
so that they both pass through the same slit of the D-slit 
aperture. Dichroic mirror DM reflects the pump beam, (b) 
The detection system. The lens IL lens is used for near field 
imaging (a;) and lens FL for measurements in the far field (p). 
BS is a 50/50% beam splitter and Dl and D2 are detectors. 



II. ENTANGLED SPATIAL QUDITS 

Entangled qudits can be produced by sending spatially 
entangled photons through D-slit apertures @, [l0| ■ In 
the near field of the slit apertures, the quantum state 
produced is well described as 0] : 

1*) = ^ E i^>ii^> 2! a) 

3 = -Id 

where \ipj) is the quantum state of a photon which passes 
through slit j and l D = (D - l)/2. For D = 2 these 
types of states have been used to identify a non-local de 
Broglie wavelength f30l ] and investigate other fundamen- 
tal aspects of non-local interference |(| H, Hf| . These non- 
local interference patterns already indicate the presence 
of entanglement in the path variables of the prepared 
two-photon states. Our aim here is to develop a formal 
technique of detecting this entanglement. 



A. Experiment 

We produced states of the form ([T]) experimentally. 
FIG. [1] shows the experimental setup, divided into two 
parts: state preparation in FIG. HJa) and the detection 
system in FIG. [IJb). FIG. HJa) illustrates the prepara- 
tion of a spatially-entangled state similar to Eq. ((T|). A 
50 mW continuous (CW) pump laser beam with 405 nm 
wavelength crosses a 2 mm Bil^Oe (BiBO) crystal and 
generates collinear photon pairs by type I SPDC. The 
photon pairs (A = 810nm) and pump beam propagate 



in the z-direction. A dichroic mirror (DM) reflects the 
pump beam and transmits the down-converted photons. 
The laser beam is focused in the center of the crystal 
by a spherical lens (SL1) with 30 cm focal length. After 
the crystal, a lens set is used to control the spatial cor- 
relation of the photons at the D-slit (D = 2,3,4) plane 
[H, [3(| , located 40 cm from the center of the crystal. A 
magnified image of the center of the crystal is projected 
at the D-slit plane by using a cylindrical lens (CL) with 
focal length /1 = 5cm and a spherical lens (SL2) with 
focal length /1 = 20cm. In this way, the correlation of 
the down-converted photons in the crystal is imaged on 
the D-slit aperture, which guarantees that the photons 
pass thr oug h the same slit, thus forming an entangled 
state [H, Hg|. In the near field just after the D-slits, the 
photons can then be described by state (JXJ) . The width 
dp of each slit is 02 = 04 = 0.08mm, 03 = 0.04mm and 
their center-to-center separation do is d-2 = d± = 0.16mm 
and d% = 0.125mm. The slits are placed perpendicularly 
to the z direction, with their smaller dimension in the 
x-direction. The length of the larger dimension of the 
slits (y-direction) is 8mm, much larger than the down- 
converted beams, and can thus be considered to be infi- 
nite. 

A detection system is set for imaging the D-slit plane 
(near field) or to project the Fourier transform of the D- 
slit (far field) at the detection plane, as shown in FIG. 
Q] (b). Both schemes use a spherical lens (IL or FL), a 
50/50% beam splitter (BS) and two single-photon detec- 
tors at the exit ports of the BS. The detectors are placed 
60 cm from the D-slit plane, and equipped with inter- 
ference filters centered at 810 nm (lOnm FWHM band- 
width) and 0.2mm diameter pinholes. For the near field 
measurement, lens IL has a focal length = 10 cm and is 
placed 12.75cm from the D-slits, giving an image magni- 
fication of M = 3.6, which was confirmed by independent 
measurement. For the far field detection, a spherical lens 
(FL) with a focal length = 30 cm in the /-/ configu- 
ration is used. We toggle between measurement schemes 
only by switching lenses IL and FL. Each detector is 
mounted in a translation stage and can be scanned in 
the x direction. Coincidences between the detectors are 
obtained with a homemade coincidence detection system 
with 5.4ns time window. 

Two-dimensional arrays of coincidence counts were ob- 
tained by scanning both detectors in the ^-direction of 
either the near- field plane or the far- field plane. FIG. [2] 
shows the coincidence distributions for D = 2,3,4 slits 
obtained in 2 seconds (near field) or 30 seconds (far field) 
for each point. Each coincidence array has more than 
20 x 20 detection points. For each near and far-field 
measurement, the detection position p x and p v in the 
x-direction of each detector was recorded, respectively. 
The position variables are then given by x = p x /M, 
and the momentum variables by p = p v j fX. The corre- 
sponding two-dimensional probability distributions were 
obtained by normalizing the coincidence distributions. 
These probability distributions will be used to evaluate 
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FIG. 2: (Color online.) Coincidence maps for (a) D — 2, 
(b) D — 3 and (c) D = 4 slits, pj (j = 1,2) is the position 
of detector Dj. Graphics on the left correspond to near- field 
(x) measurements, showing correlation at the slit planes: the 
photons pass through the same slit. Graphics on the right are 
the far-field (p) measurements, showing fringes from non-local 
interference between the two-photon wavepackets. 



the entanglement and EPR-Stccring criteria using mod- 
ular position and momentum variables presented in the 
following sections. 



III. MODULAR VARIABLES 

Let us first define the modular variables for our two- 
photon system. We will then present the entanglement 
criteria proposed in Ref. We consider position and 
momentum operators that obey the commutation rela- 
tion [x, p] = i. Choosing a scale factor £ with dimension 
of length, one can define modular variables for the con- 
tinuous x and p variables [|[ : 



x = n£ + r, 



P = m-^ + s, 



(2a) 
(2b) 



where n is the integer component of xjl and r = (x+l/2) 
mod 1-1/2 is defined so that -1/2 < r < 1/2. Similarly, 
m is the integer component of pi and s = (p + l/2£) 
mod l/£—l/2£. For the two-photon system we can define 
the global modular variables 



N± = 7ii ± n 2 , 
AI± = mi ± m,2, 



for the integer components and 



R± 
s± 



ri ± r 2 , 



(3a) 
(3b) 



(4a) 
(4b) 



for the modular remainders. 

The variances of the distributions describing the mod- 
ular variables (j = 1,2) satisfy the uncertainty relation 
| 



(A 2 r 



£ 2 (A 2 



> C = 0.078235 . 



(5) 



where the constant C is obtained numerically by calcu- 
lating the smallest eigenvalue of the operator n 2 



£H). 



Using inequality ([5]) , one can show that for any bipartite 
separable state g = J^iPiSu ® Q2i, the following inequal- 
ity holds 0: 



(A 2 N ± ) e + £ 2 (A 2 S T ) e >2C. 



(G) 



Thus, violation of this inequality indicates that the state 
is entangled. For state ((T|), one can calculate (A 2 7V_)^ = 
and £ 2 (A 2 S + )* = £ 2 /6 - (£ 2 /^)Ef=i(D - j)/Dj 2 , 
indicating a theoretical violation of the entanglement wit- 
ness © for D > 2 0. 



A. Experimental test of Variance Criteria 



Using the experimental data presented in section III Al 
and shown in Figure [TJ we tested the entanglement cri- 
teria ©. The scale factor £ is identified as the slit sep- 
aration: £ = d,D- The variances in N- and S+ used in 
the entanglement witness (|6]) were calculated using these 
probability distributions and the definitions of the modu- 
lar variables (|2|) , (|3|) and §4§ . The violation of the criteria 
(defined as the left-hand side subtracted from the right 
hand side of inequality © ) is represented by the red bars 
in FIG. [3] a). In all results presented the error bars corre- 
spond to one standard deviation (SD), obtained by error 
propagation of the Poissonian count statistics. Signifi- 
cant violation (> 3 SD) of entanglement criteria ([6]) is 
obtained only for D = 2 slits. The shaded red bars show 
the theoretical prediction for the entangled state ([]} . For 
D = 3, 4, the count rates were lower than for D = 2 slits, 
as can be noticed from the error bars, and the noise (less 
than 5%) in the near-field measurements prohibited reli- 
able identification of entanglement with the variance cri- 
teria ©, and is responsible for the discrepancy between 
the experimental values and theory. The far-field vari- 
ance (A 2 ^)^ was close to the theoretical prediction in 
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FIG. 3: (Color online.) Violation of entanglement and EPR- 
steering criteria for D — 2, 3, 4 slits. Violation is denned as 
the right-hand side of the inequality subtracted from the left- 
hand side, so that negative values correspond to violation of 
the criteria and identification of quantum correlations. The 
dark red bars show the experimental values, and the light red 
bars are the theoretical predictions calculated from the ideal 
state ([T]). Figure a) shows the variance entanglement criteria 
([6]), figure b) the entropic entanglement witness (|13[) , c) the 
variance EPR-steering criteria (|21|l and d) the entropic EPR- 
steering criteria (I24|) . For both entanglement identification 
and EPR-steering, the entropic criteria are more successful 
with our data. 



all cases. The mediocre success of the variance criteria 
for our data motivates the search for a more sensitive 
entanglement witnesses. In the next section, we present 
such a criteria. 



IV. ENTROPIC ENTANGLEMENT WITNESS 
FOR MODULAR VARIABLES 

It is known that the variance is not the optimal es- 
timate of uncertainty |37[- This has led to the devel- 
opment of entropic uncertainty relations [HI, as well as 
witnesses for entanglement (3J Htl H(| that are based 
on entropy functions. These witnesses typically outper- 
form variance-based witnesses in the case of continuous 
variable non-Gaussian states. Here we will derive an en- 
tropic entanglement witness for modular variables. Our 
derivation will follow that of Refs. [3ll |40|. We will first 
consider pure separable states, and then extend our re- 
sults to mixed states. Consider a separable pure state 



$sep(pi,P2) = $l(pi)$2(P2), 



(7) 



from which one can obtain the discrete probability distri- 
butions P(ni) and P(n2), as well as the probability den- 
sities p{s\) and p{s2) for the modular variables defined 
in Eqs. In appendix lAl we show that the modular 

variables obey the entropic uncertainty relation 



H(rij) + h(sj) > In 



1 



(8) 



where H{n) and h(s) are the discrete and continuous 
Shannon entropies [37| . From state (0 , we can calculate 
the probability densities p{S±) for global variables S± 
defined in Eq. (Q| as 



P±{Sd 



(±) 

Pi *P 2 i 



(9) 



where Pj = Pj(sj), Pj = pj(±Sj), and "*" denotes 
the convolution operation. A continuous variable whose 
probability density is the convolution of two probability 
densities satisfies the entropy power inequality [37j . and 
we can write 



h(S ± ) > iln{cxp[2/i(>i)] 



cxp[2/i(s 2 )]}. 



(10) 



To derive an entanglement witness following Refs. [3ll . 
|40| . we need a similar inequality for the discrete vari- 
ables N±. Let us assume that variables rij are non-zero 
for some finite set of values, with corresponding discrete 
probability distributions Pj . Since the quantum state ([7]) 
is pure and separable, variables n\ and n 2 are indepen- 



dent. Then, one can show that H(N±) > H(rij) [3 
That is, the sum or difference of two independent vari- 
ables increases uncertainty. It then follows that 



H(N ± ) >^ln 



1 



exp[2ff(rai)] + -exp[2H(n 2 )] 



(11) 

We note again that inequality (|TT|) is valid only for in- 
dependent variables ri\,iii. Adding inequalities (|10j) and 
l[TT|l gives 



H(N ± ) + h(S ± )>UnU ' 



2H(m)+2h{sj) 



(12) 



Using the uncertainty relation (|5J), it is straightforward 
to show that 



H{N±) +h(S±)> In ^. 



(13) 



We note that p^|) actually describes four inequalities, all 
of which are satisfied by pure separable states. 

To extend these inequalities to mixed states, we note 
that any bipartite separable state can be written as a 
convex sum of bipartite separable pure states 



(14) 



where A„, > and X w = 1. From the concavity of 
the Shannon entropy and inequality (|13[) it then follows 
that 

H(N±) e + h(S±) e >'22K > {H(N±)* w +h(S±)* w }, 

U) 

>ln^. (15) 
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Thus, inequalities (|T3)) arc satisfied by any bipartite sep- 
arable state. However, two of these are satisfied by 
all bipartite states. This follows from the fact that 
operators corresponding to N + (JV_) and S+ (S_) do 
not commute. However, it is possible to find entan- 
gled states for which H(N+) + h(S-) < ln(V2/i) or 
H(N-) + h(S+) < ln(V2/£). We thus have an entropic 
entanglement witness for modular variables. As an ex- 
ample, consider state ([T]), for which one can calculate 
H(N-) = for all D. For D = 2,3,4, we calculate 
h(S+) - ln(y/2/£) = -0.28,-0.61,-0.86, respectively, 
showing increasing violation of criteria (jT3")l with D. Wc 
note that one is free to choose the scale factor I. In ap- 
pendix IA 11 we perform a numerical investigation of I, 
and show that i — d (the slit separation) is the optimal 
choice. 



A. Experimental Test of Entropic Entanglement 
Witness 

Using the probability distributions associated to co- 
incidence counts in FIG. [2J we calculated the entropic 
modular entanglement witness (|13p for iV_ and S+ vari- 
ables. Violation of the criteria is shown by the solid red 
bars in FIG. [3]b). We obtain a significant (> 3SD) vio- 
lation for all values of D, confirming that entanglement 
is present in all cases. The shaded red bars show the 
theoretical prediction for the entangled state ((TJ). Wc 
note a considerable discrepancy between experiment and 
theory, particularly for D = 3, 4. This is due primarily 
to the noise in the near-field measurements, and also to 
the fact that the probability distribution for N_ is not 
constant, as it is for the ideal state ([1]). Nonetheless, the 
entropic entanglement witness (|13[) successfully identifies 
entanglement. 



V. EPR-STEERING CRITERIA FOR 
MODULAR VARIABLES 

One can also ask whether the correlations present in 
the state ([1]) and/or the experimental data present any 
stronger form of quantum correlation, such as EPR- 
steering correlations. 

In their seminal paper, Einstein-Podolsky and Rosen 
(EPR) discussed the appearance of a "paradox", dis- 
played by the correlations present in some bipartite quan- 
tum states. More recently, M. Rcid and collaborators 
[4TM43I ] derived a number of EPR criteria that identify sit- 
uations in which quantum correlations are strong enough 
for the EPR argument to be valid. It has been shown 
that these EPR correlations are in fact e quiv alent to 
Schrodinger's "steering" of quantum states (44(. EPR- 
steering correlations were formalized recently by Wise- 
man et al. [32| and others [H EH- It has also been 
shown that steering is related to security in certain quan- 
tum cryptography protocols [34j . Non-steerable bipartite 



states can be described by local hidden state (LHS) mod- 
els of the form [II El El : 

P{a,b) = Y,PxPM\)P PQ {b\X). (16) 

A 

Here a (b) is the outcome of measurement a (j3) and A 
is a ("hidden") variable that labels the local state, the 
distribution of which is described by classical probability 
distribution pa- -Pa(os|A) andP^Q(&|A) are the probability 
distributions for outcomes a and b for each A, where the 
subscript "Q" indicates that the conditional probability 
P/3Q(b\X) results from measurements on a quantum state. 

A. EPR-Steering Criteria with variances 

Following the derivation in Refs. [HI, E3], we consider 
the sum of variances inferred about variables of system 1 
via measurements on system 2: 

A inf (a) 2 + A inf (a') 2 (17) 

where a (a') are the outcomes of measurement a (a'), 
the inferred variance is 

A inf (a) 2 = J dbP(b) J daP(a\b)A(a\b) 2 

= J dadbP(a,b)[a- a cst ] 2 , (18) 

A(a|6) 2 is the variance of the conditional distribution 
P(a\b), and a os t is the estimate of a given outcome b 
on system 1. Let us assume that the correlations are 
described by the LHS model (fl~6|) . Then the inferred 
variance is 

A in f(a) 2 = ^p A J daP(a\X)[a - a cst ] 2 

^5> A A(a|A) 2 , (19) 

where we used / P(b\\)db = 1. Returning to expression 
(fTT]) . we have 

A inf (a) 2 + A inf (a') 2 = 5> A [A(a|A) 2 + A(a'|A) 2 ]. (20) 

Each A in the LHS model (|T6"|) labels a quantum state. 
Thus, for each A, an uncertainty principle must be satis- 
fied. Considering variables rij and Sj (j = 1,2), and us- 
ing the uncertainty relation , one arrives at the EPR- 
criteria 

A inf (n 1 ) 2 +£ 2 A inf (s 1 ) 2 > C. (21) 

Violation of criterion (|2~T|) thus indicates that the state is 
"EPR-stccrable" , and cannot be represented by the LHS 
model (|16[) . The inferred variances can be determined 
from the conditional variance using Eq. (|18[) , 
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B. EPR-Steering Criteria with Shannon entropy 

Using the LHS model (fT6|) , it is possible to show that 
this probability distribution leads to [TTj 

H(a\b) >Y,PxH Q (a\X), (22) 

A 

where the conditional Shannon entropy is H(a\b) = 
— ^2,jV{bj)H{a\b = bj). Let us consider a second set 
of conjugate observables a' and /3', for which we arrive 
at an equivalent expression. Adding the two together 
results in 

H(a\b) + H(a'\b') > ]T PA {H Q (a\X) + H Q (a'\X)} . (23) 

A 

Each A in the LHS model (fT6|) represents a different phys- 
ical realization of a quantum state. Thus, for each A, 
the state prepared should satisfy a quantum-mechanical 
uncertainty relation, such as the entropic criterion ([5]). 
Using the normalization of p\, this leads directly to an 
entropic steering criteria: 

H(m\n 2 ) + h( Sl \ S2 )>-\n£. (24) 

Any non-steerable state can be described by the LHS 
model |H]), leading directly to inequality (j2"4"]l . Viola- 
tion of inequality (|24|) thus shows that the quantum state 
shared by Alice and Bob is steerable. 



The photons are engineered to both pass through the 
same slit of a D-slit aperture, and exhibit non-local inter- 
ference effects. By introducing entanglement and EPR- 
steering criteria based on the Shannon entropy, we are 
able to detect quantum correlations for D = 2, 3, 4 slits. 
Measurement of the modular variables is the same as the 
continuous x or p variables, and thus presents no ad- 
ditional experimental challenge. The witnesses derived 
here are valid for any initial bipartite quantum state. In 
practice, the scaling parameter I can be optimized to pro- 
vide the best violation, or in the case that the slit separa- 
tion is unknown. We note that the quantum correlation 
present in both discrete and continuous variables is sim- 
ilar to correlations observed in both angle (continuous) 
and angular momentum (discrete) (48| . 

Modular variables arc a natural way to discretize 
continuous variables, and may prove extremely useful 
for continuous variable quantum information processing. 
One example application is in quantum key distribution. 
Consider a situation analogous to the experiment pre- 
sented here. Alice and Bob could measure their respec- 
tive x and p variables. When they both measure x, they 
should sec correlations in their D-dimensional discrete 
variables rt\ and n 2 , which can be used to construct a 
shared key, giving log 2 D secret bits per sifted event. Se- 
curity can be checked by using a portion of the x results 
together with the p results to test for entanglement or 
EPR-steering [34|. We expect to find similar utility in 
other quantum information tasks. 



C. Experimental Test of EPR-Steering Criteria 

We next tested the EPR-steering inequalities. Viola- 
tion of the EPR-Steering criteria (j2"Tj) and (J2"4")l . calcu- 
lated from the distributions in FIG. [21 are shown in FIG. 
[3] c) and d), respectively. In both figures, the shaded 
red bars show the theoretical prediction for the entan- 
gled state (p}. Again, the variance criteria (|2"Tj) does not 
reliably identify EPR-steering for D > 2, due principally 
to noise in the near-field measurements. On the other 
hand, the entropic criteria (|24[) detects EPR-steering for 
all values of D. Again, there is a considerable discrep- 
ancy between experiment and theory, particularly for 
D = 3,4. This is due to the noise in the near-field mea- 
surements, and also to the fact that the probability dis- 
tribution for iV_ is not constant, as it is for the ideal state 
((T|). Nonetheless, the entropic criteria (f2~4"f successfully 
identifies EPR-steering correlations. We thus expect the 
entropic entanglement criteria (|6|) and EPR-steering cri- 
teria (|2~4"|) to be quite useful for detection of correlations 
in modular variables. 

VI. CONCLUSIONS 

We have shown that entanglement in modular com- 
ponents of position and momentum variables can be de- 
tected in the spatial degrees of freedom of photon pairs. 



Appendix A: Entropic Uncertainty Relation for 
Modular Variables 

Consider a wavefunction of the form 

7 ^ = 7/r E c„exp(m0). (AI) 

* n— — oo 

The sum in (|A1|) is a Fourier series of the function 7(0), 
with discrete coefficients c„. This is analogous to an- 
gle (9) and angular momentum (n) in quantum mechan- 
ics. The probability distribution of continuous variable 
9 is 1 7 (9) 1 2 and discrete variable n is described by a dis- 
crete set of probabilities |c„| 2 , where J2 n l c »l 2 = 1- Refe- 
[H, HU derive entropic uncertainty relations for these 
variables, namely 

H(n)+h{6) >ln(27r). (A2) 

Here H(z) is the discrete Shannon entropy 

H(z) = -Y,P*^p z , (A3) 

Z 

where p z is the probability of z, h(9) is the differential 
Shannon entropy of continuous variable 9 [37j : 

TV 

h(9) = - I P(6)]nP(e)d9, (A4) 
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where P{8) is the probability distribution of 8. Note that 
all logarithms are to base e. This uncertainty relation is 
saturated for a state with one c n — 1 and the rest zero. 

This derivation was performed for a particular wave 
function. Indeed, a probability distribution describing a 
mixed state can be expanded in the form 



(A5) 



where p(8) = (8\g\p), ^(8) = (6»| 7i ), a, > and 

Y\- a,j = 1. Since this is a convex sum, and each 7,- can 
— j J 1 — 1 1 ■ j 

be expanded in the form (|A12[) . it follows [38, 49( that the 
mixed state p also obeys the entropic uncertainty relation 

Let us show that an uncertainty relation analogous to 
(|A2[) applies to the modular variables defined in Eq. . 
Following Ref. Q , consider a wave function of the form 



(A6) 



where ip n (x) is a normalized wavefunction centered at 
nd and ^ n |c„| 2 = 1. This wave function is a super- 
position of wavepackets tp, and might describe a parti- 
cle passing through an iV-slit aperture, for example (of 
course in this case all but N of the c n are zero). The 
Fourier transform of this function gives the wavefunction 
in wavevector space 



(A7) 



where 



00 

J e a ™ P Mx)dx, (A8) 



is the Fourier transform of ij) n (x). We are assuming that 
^(x) is a periodic comb of ^-functions spaced a distance 
d apart, so can write ip n (x) = ipo(x — nd), so that 4> n (j>) = 
cxp(2iTmdp)4>o(p). Then 



*(p) = ^»o(p) 2J c„exp(2i7rndp). 



(A9) 



Changing to the modular variables defined in Eqs. (|2|), 
we have 

00 

= 0o (y + «) E c« exp(2i7rnd[m/£ + s]). 



n— — 00 



(A10) 

It is natural to choose £ = d, so that exp(2i7rn^[m/^ + 
s]) = exp(2i7rn£s). If (f>o(p) ~ constant over the range 
-l/2£ to l/2£, then ^{m/i + s) sa (f> (m/£). We have 

OO 

$(p)«0o(y) 2 c n exp(2wrn&5). (All) 



In other words, $(p) ~ {m/£) g(s), where we define 

00 

g(s) = VI c„exp(2i7rn,fo). (A12) 



n— — 00 



Let us define a variable 8 = 2tt£s, so that g(8/2w£) is 
analogous to Eq. (|A1[) . We note that from —l/2£ < s < 
1/21 we have — 7r < 8 < tt. Using the uncertainty relation 
[2l then leads directly to 



ff(n) + h(s) > In 



(A13) 



By a similar argument as that outlined above for rela- 
tion (|A2[) , uncertainty relation (|A13[) applies to all quan- 
tum states p. even those that cannot be described ex- 
plicitly by wavefunction (|A12[) . The uncertainty relation 
(|A13|) is the main result of this section. We will use 
this relation to derive tests for entanglement and EPR- 
steering correlations. 



1. Choice of scaling factor 

Numerical investigation of states of the form ((T|) indi- 
cates that the optimal choice of scaling factor is I = d, 
the slit separation. One can also consider £ ^ d. How- 
ever, in this case it is necessary to consider a modified 
uncertainty relation. This is due to the fact that for £ > d 
one has coarse graining in the integer part of the x vari- 
able, and for £ < d one has coarse graining in the integer 
part of the p variable. This is equivalent to a scenario 
in which angular momentum is measured with precision 
S n > h. Bialynicki-Birula has considered these types of 
coarse grained uncertainty relations in Ref. [49|. In this 
case, the uncertainty relation (|A2[) becomes 



H(n) + H{8) > - In 



5 n Se 
2-nh' 



(A14) 



where 8 n is in units of K. Note that both entropy func- 
tions are now discrete. Construction of an entanglement 
witness using relation (|A14|) follows the same arguments 
outlined above. Explicitly, we have 



H(N±) + H{S T ) > In 



V2 



(A15) 



where 5n ± = £/d is the precision in the N± measurements 
and Ss T the precision in Sip measurements. When £ = 
d and <5s T = 0, inequality (|A15[) reduces to (TlUj) . Fig. 
|4] shows a plot of the violation of inequality (|A15|) by 
state (TT]) as a function of £/d. All results were obtained 
numerically. The largest violations occur for £ = d. 
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FIG. 4: Plot of the violation of entropic criteria (|A15|I by 
state (1} as a function of l/d for D — 2 (red circles), 3 (blue 
squares), 4 (green diamonds). The strongest violations occur 
for e/d = 1. 
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